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1. - BCS-BEC Crossover Theory and the Physical Effects of Temperature 

The study of ultracold trapped fermionic gases is a rapidly exploding subject [1-10] which is 
defining new directions in condensed matter and atomic physics. It has also captured the attention 
of physicists who study color superconducting quark matter as well as nuclear matter [11-13]. 
Indeed, it is hard, in recent times, to find a subfield of physics which appeals this broadly to 
the research community. As we come to understand the gases, experimentalists will move to 
address fermionic atoms in optical lattices; this will provide important insight to condensed 
matter physicists as analogue systems for "solving" intractable many body problems. 

What makes these gases (and lattices) so important is their remarkable tunability and con- 
trollability. Using a Feshbach resonance, one can tune the attractive two-body interaction from 
weak to strong, and thereby make a smooth crossover from a BCS superfluid to a Bose-Einstein 
condensation (EEC) [14, 15]. Not only does this allow high transition temperatures Tc (rela- 
tive to the Fermi energy Ep) but it may also provide insights into the high temperature cuprate 
superconductivity [16-19]. Furthermore, one can tune the population of the two spin states, es- 
sentially at will, in this way, allowing exploration [20-22] of exotic polarized phases such as the 
FFLO [23-25] superfluids, of interest to condensed matter, nuclear and particle physicists. One 
will be able to tune the lattice parameters such as bandwidth, on-site repulsion, even random 
disorder, etc and thereby study the famous, and as yet, unsolved Hubbard model Hamiltonian. 

This paper will concentrate on those issues relating to the effects which are of current inter- 
est to experimentalists. In particular, we will study BCS-BEC crossover in atomic Fermi gases, 
looking at a wide range of different experiments. Our group has been principally interested in 
the effects of finite temperature [26, 17, 18] and in this review, we will discuss how temperature 
T enters into the standard crossover theory, and how temperature can be measured, and how tem- 
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perature changes the character of the gas from a superfluid at low T to an unusual (but strongly 
interacting) normal fluid at moderate T, and to ultimately a Fermi gas at high T. Temperature, of 
course, plays an important role in these Fermi gas experiments. As we shall see in this Review, 
the earliest evidence for superfluidity is generally based 1 1 -8 1 on a comparison of the behavior of 
the gas at two different temperatures, presumably, one above and one below T^. Experimentally, 
one is never strictly at T = 0, and as a result there are thermal excitations of the gas which need 
to be characterized both theoretically and experimentally. 

Before we begin with the ultracold gases, it is useful to discuss to what extent the BCS- 
BEC crossover scenario relates to the high temperature copper oxide superconductors [17, 18]. 
While there is much controversy in the field of high Tc superconductivity, there is a school of 
thought 117,27-31, 19] which argues that these systems are somewhere intermediate between 
BCS and BEC. These arguments are based on the observation that is anomalously high, and 
therefore the attractive interaction is presumably stronger than in conventional superconductors. 
These are quasi-two-dimensional (2D) materials which means that they have a tendency to ex- 
hibit "pre-formed" pairs, that is, pairing in advance of Bose condensation. This is often referred 
to as fluctuation effects, but pre-formed pairs are naturaUy associated with a more BEC-Uke sce- 
nario. Importantly these pairs can be used to explain the widespread and anomalous "pseudogap" 
effects [16, 32] which are the focus [33] of most of the current attention in the field. The exis- 
tence of "pre-formed" pairs means that a characteristic (pseudogap) energy must be supplied in 
the normal state to break the pairs and create fermionic excitations. Hence we say that there is a 
gap or, more precisely, a pseudogap in the fermionic excitation spectrum. 

The case that the crossover scenario is relevant to the cuprates was made quite eloquently 
by A. J. Leggett in a recent status report [19] on high superconductivity. In this article he 
summarized the eight salient "facts" about the cuprates. It is worth quoting one here, which 
relates to their anomalously short coherence length. "The size of the pairs is somewhere in 
the range 10-30 A — from measurements of the upper critical field, Fermi velocity and T^. This 
means that the pair size is only moderately greater than the inter-conduction electron in-plane 
spacing, putting us in the intermediate regime of the so-called Bose-Einstein condensate to BCS 
superconductor (BEC-BCS) crossover, and leading us to expect very large effects of fluctuations 
(they are indeed found)." 

The field ofBCS-BEC crossover is built around early observations by Eagles [15] and Leggett 
[14] that the BCS ground state, proposed by Bardeen, Cooper, and Schrieffer in 1957, is much 
more general than was originally thought. If one increases the strength of the attraction and self- 
consistently solves for the fermionic chemical potential p,, this wave function will correspond 
to a more BEC-like form of superfluidity. Knowing the ground state what is the nature of the 
superfluidity at finite T? That is the central question we will address in this Review. 

Even without a detailed theoretical framework we can make three important observations. 

• As we go from BCS to BEC, pairs will form above Tc without condensation. In the normal 
state, it pays in general to take advantage of the attractive pairing interaction. Only in the 
extreme BCS limit do pairs form exactly at Tc. 

• The fundamental statistical entities in these superfluids are fermions. We can think of 
pairs of fermions as a form of "boson", although the statistics are not precisely the same. 
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Figure 1. - Contrasting beliavior of tlie excitation gap A(T) and superfluid order parameter A^dT) versus 
temperature. The height of the shaded region roughly reflects the density of noncondensed pairs at each 
temperature. 



We measure these "bosonic" or pair-degrees of freedom indirectly via the fermionic gap 
parameter A(r). In the fermionic regime this parameter is the minimum energy which 
must be supplied to create fermionic excitations. It tells us about bosons indirectly through 
the binding together of two fermions. 

• In general there will be two types of excitations in these BCS-BEC crossover systems. 
Only in strict BCS theory are the excitations of solely fermionic character, and only in the 
strict EEC limit will the excitations be solely of bosonic character More generally in the 
intermediate case (often called the "unitary" regime) the excitations consist of a mix of 
both fermions and bosons. 

These observations are illustrated by Figs.[2and|2] In Fig.[2we schematically plot the gap 
parameter A(T) as a function of T, along with the superfluid order parameter Asc{T). The 
former, which represents the "bosonic" degrees of freedom, shows that pairs continuously form 
once temperature is less than a crossover temperature T*, while the order parameter turns on as 
in a second order phase transition at Tc- The height of the shaded region reflects the number 
of noncondensed pairs. This number increases monotonically with decreasing T, until Tc is 
reached. A T further decreases below Tc the number of noncondensed pairs begins to decrease 
monotonically due to the condensation of zero momentum pairs. 

In Fig. 12] we present a schematic plot of the excitation type, which shows that between BCS 
and EEC (i.e., in the unitary regime) there will be a mix of fermions and bosons. These bosons 
and fermions are not separate fluids, but rather they are strongly inter-connected. Indeed, the gap 
in the fermionic spectrum (related to A) is a measure of the number of bosons in the system. 

2. - Theory Outline 

In this paper we will presume the ground state is given by the standard BCS-like wavefunc- 
tion. 
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Figure 2. - Schematic illustration of excitations in the BCS, unitary and BEC regimes. The black discs 
represent fermionic excitations. Pair excitations become progressively dominant as the system evolves 
from the BCS to BEC regime. 



where cj^. ^ and Cj^ ^ are the creation and annihilation operators for fermions of momentum k and 
spin fT =t, i- The variational parameters Uk and Wk are usually represented by the two more 
directly accessible parameters Asc{T = 0) and /i, which characterize the fermionic system. 
Here Asc{T = 0) is the zero temperature superconducting order parameter. These fermionic 
parameters are uniquely determined in terms of the attractive interaction U and the fermionic 
density n. The variationally determined self consistency conditions are given by two BCS-like 
equations which we refer to as the "gap" and "number" equations, respectively. 

At T = 0, the effects of BEC-BCS crossover are most directly reflected in the behavior of 
the fermionic chemical potential /i, which is associated with the number equation. We plot the 
behavior of fi in Fig.|3] which indicates the BCS and BEC regimes. In the weak coupling regime 
H = Ep and ordinary BCS theory results. With increasing \U\, fj, begins to decrease, eventually 
crossing zero and then ultimately becoming negative in the BEC regime. We generally view 
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Figure 3. - Typical behavior of the chemical potential /i at T = in the three regimes, as a function of 
the interaction strength U /Uc, or, equivalently, l/kpa. As U /Uc increases from 0, the chemical potential 
/i starts to decrease and then becomes negative. The character of the system changes from fermionic (/i > 
0) to bosonic (jj, < 0). The pseudogap (PG) or unitary regime corresponds to non-Fermi liquid based 
superconductivity, and Uc{< 0) corresponds to critical coupling for forming a two fermion bound state in 
vacuum [17], i.e., the unitary point where the two-body s-wave scattering length a diverges. 
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Table I. - BCS theory by way of BEC analogy. Here we compare condensation in composite and point 
bosons; fiB is the bosonic chemical potential. No and Nt are the number of condensed and noncondensed 
bosons, respectively. We define fipair as the chemical potential for the noncondensed pairs. Here A(r) is 
the total fermionic gap which contains contributions from the noncondensed (Apg) and condensed terms 
f j- In the strict BCS limit Apg — 0, so that the order parameter and gap are identical. 



Composite bosons 



Point bosons 



Pair chemical potential 



Impair — 0, T ^ Tc 

Leads to BCS gap equation for A(T) 



Total "number" of pairs 



N ^ No + Nt 



Noncondensed pairs 



/i = as a crossing point. For positive /i the system has a remnant of a Fermi surface, and we 
say that it is "fermionic". For negative fi, the Fermi surface is gone and the material is "bosonic". 

One can debate whether other ground states ought to be considered. Indeed the work of 
the Camerino group [34, 35] is based on a finite T approach first introduced by Nozieres and 
Schmitt-Rink (NSR) [36, 31]. This leads to a different, and not as readily characterized ground 
state. We list some arguments in support of the ground state in Eq. i2.1\ . (i) This is the basis for 
the widely studied Bogoliubov-de Gennes approach, which can be applied [37-40] to the BCS- 
BEC crossover problem at T = 0. (ii) At arbitrarily strong coupling (and T — 0), this ground 
state can be shown to coincide with a Gross-Pi taevskii (GP) description [41] of the boson system, 
(iii) This ground state is the basis for the rapidly proliferating theoretical literature [40,42-47] on 
spin polarized Fermi superfluids. In addition to our own work there have been some additional 
studies which include the effects of temperature [48-50], albeit at a lower order mean field theory 
than considered here. 

We begin at the more physical level by stressing the analogy between condensation in this 
composite boson or fermionic superfluid and condensation in a gas of ideal point bosons. Our 
microscopic theory treats self-consistently two-particle and one-particle Green's functions on an 
equal footing. Because the physics is so simple and clear, we can fairly readily anticipate the 
form of the central equations of this BCS-BEC generalization of BCS theory. It is important to 
stress, however, that these equations can be derived more rigorously from a truncated series of 
equations of motion for the appropriate Green's functions [16,51]. 

There are three principal equations which govern Bose condensation: the vanishing of the 
bosonic chemical potential at all T < Tc is the first. We will refer to this condition as the "BEC 
condition". It is related to the usual Thouless criterion for superconductivity, but the latter is 
generally associated only with the temperature Tc. The second equation is the boson number 
equation. All "bosons" must be accounted for as either condensed or noncondensed. The third 
equation is the number of noncondensed "bosons", which are created by thermal excitations. 
This is determined simply by inserting the known excitation spectrum of the excited pairs or 
bosons, into the Bose distribution function b{x). With this equation, and the second equation. 
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one can then deduce the number of condensed bosons. 

These three central equations for bosons are indicated in Table I, on the far right, for true 
point bosons, and in the second column for the composite bosons which appear in fermionic 
superfluids. For these composite bosons the quantity which provides a measure of the "number" 
of bosons (N) is given by (T) (up to a constant coefficient, Z). How does one quantitatively 
establish the appropriate "boson number" for the fermionic case? This is determined via the self 
consistent gap equation for A(T), which, in turn, is determined using the first condition: that 
the pair chemical potential fipair is zero at and below T^.. How does one compute the number of 
excited pairs? Once the gap equation is interpreted in terms of the appropriate noncondensed pair 
propagator (see below), then one knows the related excitation spectrum fig of this propagator. 

The quantity Z which appears in the last equation of the Table (for the composite bosons) 
gives the relation between the gap associated with noncondensed pairs (A^^) and the number of 
all pairs [J^ K'^q)^- It can be readily calculated in this theory; once one has the noncondensed 
pair propagator, Z appears as the inverse residue. (Deep in the BEC regime [52], Z is relatively 
simple to compute, for here the boson number density approaches the asymptote n/2). More 
precisely, the total number of bosons in the present case has to be determined self-consistently 
through the gap equation for A(T). It also involves the fermion number equation through the 
related fermionic chemical potential. 

To be consistent with the ground state variational conditions, the vanishing of the pair chemi- 
cal potential is associated with a particular form for the pair propagator involving dressed Green's 
functions. These, in turn, determine the fermionic chemical potential through the fermion num- 
ber equation. 

21. Microscopic T-matrix Scheme. - Next, we implement this picture microscopically via a 
T-matrix approximation. We include spin indices throughout so that it will be clear how to apply 
this scheme to spin polarized superfluids. This means that we consider the coupled equations 
between the particles (with propagator G) and the pairs [with propagator t{P)] and drop all 
higher order terms. This theory does not include direct "boson-boson" interactions, although 
the pairs do interact indirectly via the fermions, in an averaged or mean field sense. Here, for 
all T < Tc, the BEC condition is interpreted as requiring that the pair chemical potential iJ,pair 
associated with the noncondensed pairs, vanish. Within a T-matrix scheme, the pair propagator 
is given by 



where x is the pair susceptibility. The function x(-P) is, in many ways, the most fundamental 
quantity we introduce. We will show that one obtains consistent answers between T-matrix 
based approaches and the BCS-Leggett ground state equations, provided the components of the 
pair susceptibihty, defined by 



(2.2) 



(2.3) 



x(^') = ^[xu(^)+XiT(^')] 
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are given by the product of one dressed and one bare Green's function 

(2.4) xniP) = J2'^oi{P~K)Gi{K), xil{P)=T.^oi{P-K)G^{K), 

K K 

where P = {ifti, p), and G and Go are the full and bare Green's functions respectively. Here 
Gq]^{K) — iujn — ^k,o-, ^k.cr = — l^a, ^k = h^k^ /2m is the kinetic energy of fermions, and 
ficr is the fermionic chemical potential for spin cr =t, i. Throughout this paper, we take h ^ 1, 
ks — 1, and use the four-vector notation K = (icjn, k), P = {ifli, q), J^k = P En X^k' 
where w„ = (2n + l)7rT and il/ — 21ttT are the standard odd and even Matsubara frequencies 
[53] (where n and I are integers). 

We now evaluate the BEC condition 



(2.5) <-i(0) = = C/-i + x(0). 
The one-particle Green's function for fermions with spin a is 

(2.6) G-\K) = G^^iK) - Y.,{K) = zc^„ - Cfc<x - ^a{K) 
where a = —a, and the self-energy is of the BCS-like form 



(2.7) 



It should be noted that we use a contact potential so that the symmetry factor (p\^ associated with 
the pairing interaction is trivially (^j. = 1. For a nontrivial </Jk, one only needs to replace A 
with A(^k in Eq. ( 12. 7> . We will see below how this form for the self energy very naturally arises 
(below Tc) in a T-matrix approach. Thus 



(2.8) 



A^ 



(2.9) 



Now we are in position to calculate the pair susceptibility at P = 

1 



x(0)=XTl(0)=XiT(0) = -5] 



Substituting this expression into our BEC condition Eq. i2.5\ . we obtain the gap equation 



(2.10) 



I- f{Eui)- f{Ek^) 
2Ek 



1 - 2f{Ek) 
2Ek 



after carrying out the Matsubara summation. Here /i = + i-i i)/2 and h = — Mi)/2, 
Ek = ■\/^^ + A^, Ek^ ~ —h + Ek and Eki = h + Ek, where — — fi. In addition, 
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we define the average f{x) = [f(x + h) + f{x — h)]/2, where f(x) is the Fermi distribution 
function. 

The coupUng constant U can be replaced in favor of the dimensionless parameter, 1/kpa, 
via the relationship m/ (Ana) = 1/U + J^ki'^^k)^^ ^ where a is the two-body s-wave scattering 
length, and fci? is the noninteracting Fermi wave vector for the same total number density in the 
absence of population imbalance. Therefore the gap equation can be rewritten as 



Here the "unitary scattering" limit corresponds to resonant scattering where a —>■ oo. This 
scattering length is tunable by magnetic field application and we say that we are on the BCS or 
BEC side of resonance, depending on whether the fields are higher or lower than the resonant 
field, or alternatively whether a is negative or positive, respectively. 

Finally, in terms of Green's functions, we readily arrive at the number equations: = 
J2k Ga{K), which are consistent with their ground state counterparts 



where the coherence factors = (1 ± ^k/£'k)/2. 

Next we use this T-matrix scheme to derive Eq. M.ll and separate the contribution from 
condensed and noncondensed pairs. The diagrammatic representation of our T-matrix scheme 
is shown in Fig.|4] The first line indicates tpg, and the second the total self energy. One can see 
throughout the combination of one dressed and one bare Green's function, as represented by the 
thick and thin lines. The self energy consists of two contributions from the noncondensed pairs 
or pseudogap (pg) and from the condensate (sc). There are, analogously, two contributions in 
the full T-matrix 



(2.11) 




(2.12) 




k 



(2.13) 



t 



(2.14) 



U 



(2.15) 



tsc{P) 



1 + UxiP)' 



where we write A^c = I]k(c-kiCkT)- 

Similarly, we have for the fermion self energy 



(2.16) J:^{K) = K'{K) + ^Pa{K) = ^ t{P)Go.AP - K) . 



p 



We can see at once that 



(2.17) Y.%-{K) = tsc{P)GoAP -K)^ -GoA^K)Al 

p 
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The vanishing of the pair chemical potential implies that 

(2.18) t;/(0) = + x(0) = 0, r < . 

Moreover, a vanishing chemical potential means that tpg (P) diverges P = 0. Thus, we may 
approximate [54] Eq. ( I2.16> to yield 

(2.19) I],(i^)«-Go,ff(-i^)A2, 
where 

(2.20) A'{T)^Al{T) + Alg{T), 
Importantly, we are led to identify the quantity Apg 

(2.21) Ks^-T. ^PaiP)- 

Note that in the normal state (where impair is non-zero) Eq. (I2.19> is no longer a good approxi- 
mation. 

We now have a closed set of equations for addressing the ordered phase. This approach can be 
readily generalized [55] to treat more exotic polarized phases such as the LOFF state [23,25,24]. 
We can similarly extend this approach to temperatures somewhat above Tc, by self consistently 
including a non-zero pair chemical potential. This is a necessary step in addressing a trap as 
well [56]. Additionally, the propagator for noncondensed pairs can now be quantified, using 
the self consistently determined pair susceptibility. At small four-vector P, we may expand the 
inverse of tpg, after analytical continuation (iili + iO+), to obtain 

(2.22) tpg^ ^ain^ + z(^n- ^ + fipa.r + iTp^ , 

where the imaginary part Tp rapidly as p ^ below Tc- Because we are interested in the 
moderate and strong couphng cases, we drop the aifl^ term in Eq. ( I2.22K and hence 

(2.23) tpg (P) = „ „ / ' ^ .„ , 

l L -\- j-lpair i ^-L p 

where we associate 

(2.24) Qp- ^ 



2M* 

This establishes a quadratic pair dispersion and defines the effective pair mass, M*. They can be 
calculated via a small p expansion of x(P), 

dx 1 1 d'^x 
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Figure 4. - T-matrix and self-energy diagrams for tire present T-matrix sclieme. Tiie self-energy comes 
from contributions of both condensed (Esc) and noncondensed (Spg) pairs. Note that there is one dressed 
and full Green's function in the T-matrix. Here tpg represents the propagator for the noncondensed pairs. 



Finally, one can rewrite Eq. ( 12. 2U as 

(2.26) A2^(T) = Z-i^6(f7p). 

p 

We now return to the strong analogies between our T-matrix approach and Bose condensation 
of point bosons, as summarized in Table I. We have three central equations. 

1 . The pair chemical potential must vanish at and below Tc, 

(2.27) Mpa.r =0, (T < T,). 

Importantly this condition leads to the mean field gap equation derived in Eq. ( I2.10> . 

2. There must be a conservation of the total number of (composite) "bosons" in the system. 
For this condition, our central equation is Eq. ( I2.20> . Here it is understood that the number of 
"bosons" is effectively represented by the parameter A^(T). Unlike the point boson case, the 
"total boson number" is temperature dependent and has to be self-consistently determined. 

3. The number of noncondensed pairs is readily computed in terms of the pair dispersion, 
just as in conventional BEC. For this condition our central equation is Eq. ( I2.26> . Then, just as 
in conventional BEC, the number of condensed bosons (proportional to A^^) is determined by 
the difference between A^(T) and Apg{T). This, in turn, determines Tc as the lowest temper- 
ature(s) in the normal state at which noncondensed pairs exhaust the total weight of A^ so that 
Apg = A^. Solving for the "transition temperature" in the absence of pseudogap effects [48-50] 
leads to the quantity T*^^. More precisely, T^^^ should be thought of as the temperature at 
which the excitation gap A(T) vanishes. This provides a reasonable estimate, for the pairing 
onset temperature T*, (when a stable superfluid phase exists). This is distinguished from the 
transition temperature. We note that T* represents a smooth crossover rather than an abrupt 
phase transition. 

It should be stressed that the dispersion relation for the noncondensed pairs is quadratic. 
While one will always find a linear dispersion in the collective mode spectrum [57], within the 
present class of BCS-BEC crossover theories, the restriction to a T-matrix scheme means that 
there is no feedback from the collective modes onto the pair excitation spectrum. In effect, the 
T-matrix approximation does not incorporate pair-pair interactions at a level needed to arrive 
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1//cpa ^/kpa 

Figure 5. - Typical beiiavior of Tc as a function of l/kpa in a homogeneous system (left panel) and 
in a trapped Fermi gas (right panel). Tc follows the BCS predictions and approaches the BEC asymptote 
0.218rF and 0.518Tf in the homogeneous and trapped cases, respectively. In contrast to the homogeneous 
case, the BEC asymptote in a trap is much higher due to a compressed profile for trapped bosons. In the 
homogeneous case, Tc reaches a maximum around l/kpa — and a minimum around where fi = 0. In the 
trapped case, this maximum/minimum behavior is washed out largely by the shrinking cloud size as l/kpa 
increases. 



at this expected linear dispersion in the pak excitation spectrum. Nevertheless, this level of 
approximation is consistent with the underlying ground state wavefunction. 

3. - Behavior of Tc and Trap Effects 

Before turning to experiment, it is important to discuss the behavior of Tc which is plotted as 
a function of scattering length in the left panel of Fig.|5lfor the homogeneous case, presuming 
s-wave pairing. Starting from the BCS regime this figure shows that Tc initially increases as the 
interaction strength increases. However, this increase competes with the opening of a pseudogap 
or excitation gap A(Tc). Technically, the pairs become effectively heavier before they form 
true bound states. Eventually Tc reaches a maximum (very near unitarity) and then decreases 
slightly until field strengths corresponding to the point where /i becomes zero. At this field value 
(essentially where Tc is minimum), the system becomes a "bosonic" superfluid, and beyond this 
point Tc increases slightly to reach the asymptote 0.218Ep corresponding to an ideal Bose gas. 

Trap effects change these results only quantitatively as seen in the right panel of Fig.|5] Here 
the maximum in Tc may no longer be visible. The calculated value of Tc (~ 0.3Tp) at unitarity 
is in good agreement with experiment [8, 58] and other theoretical estimates [59]. To treat these 
trap effects one introduces the local density approximation (LDA) in which Tc is computed under 
the presumption that the chemical potential jj. ^ fi — V{r) . Here we consider a spherical trap 
with V{r) ~ ^muj'^r^. The Fermi energy Ep is determined by the total atom number N via 
Ep = ksTp = huj{3Ny^'^ = h'^kp/2m, where kp is the Fermi wavevector at the center of 
the trap. It can be seen that the homogeneous curve is effectively multiplied by an "envelope" 
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Figure 6. - Typical spatial profile for density n(r) and fermionic excitation gap A{r) of a Fermi gas in 
a trap. The curves are computed at T = Tp/i and at unitarity, where l/kpa — 0. Here Rtf is the 
Thomas-Fermi radius, which gives the cloud size at T = in the noninteracting limit. 



curve when a trap is present. This envelope, with a higher BEC asymptote, reflects the fact that 
the particle density at the center of the trap is higher in the bosonic, relative to the fermionic 
case. In this way Tc is relatively higher in the BEC regime, as compared to its counterpart in the 
homogeneous case. 

Figure|6]is a central one, for it prepares us for understanding various experiments. It presents 
a plot of the position dependent excitation gap A(r) and particle density n{r) profile over the 
extent of the trap. An important point needs to be made: because the gap is largest at the center 
of the trap, bosonic excitations will be dominant there. At the edge of the trap, by contrast, 
where fermions are only weakly bound (since A(7-) is small), the excitations will be primarily 
fermionic. We will see the implications of these observations as we examine thermodynamic [60] 
and radio frequency (RF) spectra data [61] in the ultracold gases. 

One should appreciate that temperature is not straightforward to measure in these cold gases; 
it is most readily measured at the Fermi gas endpoint [3] or in the deep BEC regime [62]. At 
unitarity, the physical temperature can be extracted using phenomenological fits to the particle 
density profiles based on the universality hypothesis [63,58,64] with proper re-calibration [8,65] 
below Tc- For a more general magnetic field one has to resort to adiabatic sweep thermometry 
[60,66]. Here, the magnetic field of interest is accessed via a slow, adiabatic, or isentropic, 
magnetic field sweep starting from either the BCS or BEC endpoints, where the temperature 
(and the entropy [67]) are known. A finite temperature theory of BCS-BEC crossover is required 
to calculate the entropy [60] at general magnetic fields. In this way, the physical temperature 
can be associated with the measured endpoint temperature. Indeed, we will see below that the 
temperature which appears in the measured superfluid phase diagram [3] or in the RF pairing gap 
experiments [62] is given in terms of the endpoint temperature. 
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Figure 7. ~ Decomposition of density profiles at various temperatures at unitarity. Here green (light gray) 
refers to the condensate, red (dark gray) to the noncondensed pairs and blue (black) to the excited fermionic 
states. Tc — 0.27Tf, and Rtf is the Thomas-Fermi radius. The presence of noncondensed pairs is 
essential [65] for explaining why there are no sharp features in these profiles, associated with the interface 
of the normal and superfluid regions. Here Us, Up, and ngp denote density of superfluid, incoherent pairs, 
and fermioinic quasiparticle, respectively. 



4. - Experimental Evidence for a Pseudogap in Cold Gases 

Our finite temperature generalization of tlie BCS-like ground state has introduced the concept 
of a "pseudogap". This pseudogap in the fermionic spectrum should be viewed as synonymous 
with the concept of noncondensed pairs, or with pairs which have a finite center of mass mo- 
mentum. They are important both above and below Tc- In this section we want to explore the 
evidence for these noncondensed pairs using three different experiments: density profiles, nor- 
mal state thermodynamics and RF pairing gap spectroscopy. In this section we will consider the 
case of an unpolarized gas. 

In Fig. 0we plot a decomposition of the particle density profiles [65] for various temperatures 
above and below Tc- The various color codes (or gray scales) indicate the condensate along with 
the noncondensed pairs and the fermions. This decomposition is based on the superfluid density 
so that all atoms participate in the condensation at T 0. 

An important observation should be noted. The noncondensed pairs are responsible for 
smoothing out what otherwise would be a discontinuity between the fermionic and condensate 
contributions. This leads to a featureless profile, in agreement with experiment [68,69]. Indeed, 
these experimental observations presented a challenge for previous theories [70, 63] which ig- 
nored noncondensed pairs, and therefore predicted an effectively bimodal profile with a kink at 
the edge of the superfluid core. One can see from the figure that even at Tc, the system is different 
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Figure 8. - Pseudogap effects as apparent from thermodynamics. From Ref. [8]. The fact that the exper- 
imental data (symbols) (and the two theoretical curves) for noninteracting and unitary Fermi gases do not 
merge until higher T* > Tc is consistent with the presence of a normal state pseudogap. 



from a Fermi gas. That is, noncondensed pairs are present in the central region of the trap when 
the condensate is gone. Even at T /Tc = 1.5 there is a considerable fraction of noncondensed 
pairs. It is not until around T* = 2Tc for this unitary case, that noncondensed pairs have finally 
disappeared. 

We next turn to a detailed comparison of theory and experiment for thermodynamics. Figure|8l 
presents a plot of energy i? as a function of T comparing the unitary and non-interacting regimes. 
The solid curves are theoretical while the data points are measured in ^Li [8]. There has been a 
recalibration of the experimental temperature scale in order to plot theory and experiment in the 
same figure. The latter was determined via Thomas-Fermi fits to the density profiles. To arrive at 
the caUbration, we applied the same fits to the theoretically produced density profiles, examples 
of which appear in Fig.0 Good agreement between theory and experiment is apparent in Fig.|8] 
In the figure, the temperature dependence of E reflects primarily fermionic excitations at the 
edge of the trap, although there is a small bosonic contribution as well. Importantly one can see 
the effect of a pseudogap in the unitary case. The temperature T* is visible from the plots as that 
at which the non-interacting and unitary curves merge. This corresponds roughly to T* « 2Tc. 

Measurements [62] of the excitation gap A can be made more directly, and, in this way 
one can further probe the existence of a pseudogap. This pairing gap spectroscopy is based on 
using a third atomic level, called |3), which does not participate in the superfluid pairing. Under 
application of RF fields, one component of the Cooper pairs, called |2), is excited to state |3). 
If there is no gap A then the energy it takes to excite |2) to |3) is the atomic level splitting 1^23 ■ 
In the presence of pairing (either above or below Tc) an extra energy A must be input to excite 
the state |2), as a result of the breaking of the pairs. Figure |9] shows a plot of the spectra for 
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Figure 9. - Experimental RF Spectra for "Li at unitarity at 834 G. The temperatures labeled in the figure 
were computed theoretically at unitarity based on adiabatic sweeps from BEC. The two top curves, thus, 
correspond to the normal phase, thereby, indicating pseudogap effects. Here Ef = 2.5/iK, or 52 kHz. 
FromRef. [62]. 

^Li near unitarity for four different temperatures, which we discuss in more detail below. In 
general for this case, as well as for the BCS and BEC limits, there are two peak structures which 
appear in the data and in the theory [71,61]: the sharp peak at 1^23 = which is associated 
with "free" fermions at the trap edge and the broader peak which reflects the presence of paired 
atoms; more precisely, this broad peak derives from the distribution of A in the trap. At high 
T (compared to A), only the sharp feature is present, whereas at low T only the broad feature 
remains. The sharpness of the free atom peak can be understood as coming from a large phase 
space contribution associated with the 2^3 excitations [61]. These data alone do not directly 
indicate the presence of superfluidity, but rather they provide strong evidence for pairing. 

It is interesting to return to discuss the temperatures in the various panels. What is measured 
experimentally are temperatures T' which correspond to the temperature at the start of an adia- 
batic sweep from the BEC limit to unitarity. Here fits to the BEC-like profiles are used to deduce 
T' from the shape of the Gaussian tails in the trap. Based on knowledge [60] about thermo- 
dynamics (energy E in the previous figure or, equivalently, entropy S), and given T', one can 
then compute the final temperature in the unitary regime, assuming S is constant. Indeed, this 
adiabaticity has been confirmed experimentally in related work [69]. We find that the four tem- 
peratures are as indicated in the figures. Importantly, one can conclude that the first two cases 
coiTespond to a normal state, albeit not far above Tc. In this way, these figures suggest that a 
pseudogap is present as reflected by the broad shoulder above the narrow free atom peak. 

5. - Establishing Superfluidity in Cold Fermi Gases 

From the time of the earliest discoveries [72, 1-8] there was a drive to establish the existence 
of superfluidity which is more difficult on the BCS than on the BEC side [2, 1] of the resonance. 
There has been a sequence of experiments which have effectively made this case, beginning 
first with fast sweep experiments [3,4], then thermodynamical measurements [8], and finally 
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Figure 10. - (Color) Earliest evidence for superfluidity: phase diagram of^K as a function of (T/Tf) 
and 1/kpa. This compares both experiment (left panel) and theoretical computations (right panel). In the 
left panel, a contour plot of the measured condensate fraction Nq/N as a function of 1/kpa and effec- 
tive temperature {T /Tp)'^ is compared with theoretically calculated contour lines at Ns/N — (at the 
superfluid transition, black curve) and 0.01 (red curve). The experimental data have an overall systematic 
uncertainty of approximately 0.1 in 1/k^a. The experimental contour at No/N — 0.01 and the theo- 
retical line for Ns/N — 0.01 are in good agreement. The dashed line represents the naive BCS result 
Tc/T]^ « O.eiSe'"/^*^". The right panel represents a more complete theoretically computed equilibrium 
phase diagram, with contour lines for Na/N. Here all temperatures are measured in the Fermi gas regime. 
FromRef [73]. 



detection of quantized vortices [5]. We discuss the first two methodologies here in the context of 
our theoretical framework. We limit our discussion first to the case of unpolarized superfluids. 

The left panel in Fig.^Jis a plot of the first phase diagram representing the condensate frac- 
tion vs 1 /kpa, as obtained in Ref. [3] for ^"K. Subsequently, similar studies [4] were undertaken 
for ^Li. The figure is based on starting the system off in the free Fermi gas regime where it can be 
associated with an initial known temperature {T/Tp)^ , and then adiabatically sweeping closer 
to unitarity. [Here Tp = Tp, kp = kp and [T/Tp]^ are all measured in the noninteracting 
Fermi gas limit]. Once the near-unitary gas (of fixed entropy) is obtained, a fast sweep is made 
to the BEC regime, where the condensate fraction can be read off from a bi-modal profile. The 
presumption here, for which there is considerable experimental support [3,4], is that, even if the 
condensate fraction is not conserved upon a fast sweep to BEC, the presence or absence of a 
condensate will be preserved. The time frame for the sweep will not allow a condensate to form 
in the BEC if there were none present near unitarity, nor will it allow a condensate to disappear 
if it was present initially. 

The lines drawn on top of the experimental contour plots are the calculated [73] condensate 
fraction contours as a function of the adiabatic sweep-projected temperatures [T/Tp)'^ for a 
0% and 1 % condensate fraction. These essentially correspond to the normal-superfluid phase 
boundary which is expected to be rather well measured in these sweep experiments. The figure to 
the right presents a full plot of the theoretical condensation fraction, importantly, measured with 
respect to the adiabatic sweep-projected temperatures {T/Tp)^ . This, then, is the theoretical 
phase diagram. For the 1% case, the overall trends yield good agreement between theory and 
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Figure 11.- Evidence for a superfluid phase via plots of the energy E vs physical temperature T. The 
upper curve (long dashed line) and data points correspond to the BCS or essentially free Fermi gas case, 
and the lower curve and data correspond to unitarity. The latter provide indications for a phase transition 
via a slope change. The thin dashed line represents a power law fit to the unitary data below the transition. 
The inset shows how temperature must be recalibrated below Tc. Here /3 is a parameter which characterizes 
the interaction energy at unitarity for T = 0, and T is a empirical temperature scale. From Ref. [8]. 



experiment, except for the small "overshoot" (of unknown origin) which appears in the BEC 
side of the data. 

The second generation experiments which helped to establish superfluidity were based on 
thermodynamical measurements. In Fig. ^2 we show data from Ref. [8]. What is plotted is the 
measured energy as a function of temperature on a log-log scale. This temperature represents 
a theoretical recalibration of an effective profile-measured temperature T . The key feature here 
is that the data (indicated by the points) show an abrupt change at a temperature one can call 
Tc. This abrupt change occurs for the unitary scattering case. No such feature is seen for the 
noninteracting Fermi gas, also plotted in the figure. Furthermore, this abrupt change in E vs 
T is seen even in the raw data points (not shown here), without appealing to a theoretical re- 
calibration. All of this is very suggestive of a specific heat discontinuity, which is to be associated 
with a phase change - presumably to a superfluid phase. 

The inset in the figure shows how the effective temperature T which is obtained from a 
Thomas-Fermi fit to the calculated density profile compares with the physical temperature T. 
The inset shows that a recalibration is necessary only below Tc, to account for the presence of a 
condensate. We stress that the observation of a phase change made by the Duke group [8] is not 
dependent on this recalibration. Hence these experiments provide good evidence for a transition 
between a normal and superfluid phase. 

The last generation experiment to make the case for superfluidity was the observation of 
quantized vortices by the MIT group [5]. These very convincing experiments will be discussed 
elsewhere in this lecture series. 
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6. - Fermi Gases with Imbalanced Spin Population 

The latest excitement in the field of trapped fermions pertains to gases with deliberately im- 
balanced spin populations [20-22]. In large part this is motivated by interest from theorists in 
other disciplines such as dense QCD and (isospin asymmetric) nuclear matter [11-13, 74]. From 
the condensed matter viewpoint there has been an underlying interest in exotic and intriguingly 
elusive phases associated with Zeeman effects in superconductors, such as that proposed by 
Larkin and Ovchiimikov and by Fulde and Ferrell (LOFF). In the LOFF state [23] the conden- 
sate has a net momentum of a fixed q or of -|-q and — q. Even more elaborate crystalline lattices 
of various have also been contemplated [24]. An additional, and very important motivation 
for these studies is associated with the recent claims [22] that when there is a population imbal- 
ance, the density profiles will indicate whether or not superfluidity is present and they, moreover, 
provide an internal mechanism for thermometry. Thus, because these experiments are claimed 
to identify Tc itself, and because they pertain to thermometry, it should be clear that a theoretical 
analysis of these experiments requires an understanding of the effects of T 0. 

We begin our discussion by summarizing some key experimental observations [20-22]. In a 
trapped cloud there appears to be a superfluid core, which, at the lowest temperatures, is unpo- 
larized. Outside of this core there is a normal region, in which both spin components are present 
and this carries a significant fraction of the polarization. Beyond this "mixed, normal region" 
there is a free Fermi gas composed of the majority spin species, which carries additional polar- 
ization. At suffiently low temperatures, there appears to be a form of phase separation in which 
the superfluid and normal phases are associated with zero and finite polarization, respectively. 

We turn now to a theoretical understanding [56, 75] of these experiments atT ^ first for 
the case of a homogeneous gas, and then later in the trapped configuration. There has been some 
work along these same lines elsewhere in the literature [49], although without incorporating 
noncondensed pair effects. We focus on the Sarma or breached pair phase [76]. The figures we 
present do not include the more exotic LOFF or phase separated states. The formal structure 
for addressing the former [55] is very similar to that of the Sarma state; we will defer a brief 
discussion of these until the end. To help with the clarity of the presentation, we state our major 
conclusions for the homogeneous and trapped configurations in the unitary regime, at the outset. 

• In the homogeneous case temperature serves to stabilize the polarized superfluid phase. 
Such a phase is unstable in the ground state. 

• In a trap at low T, the superfluid core remains unpolarized. Except at very low T, spin 
polarization will be mostly carried by fermions coexisting with strongly interacting, ther- 
mally excited noncondensed pairs which appear outside the core. Thus, pseudogap effects, 
which have been emphasized throughout, are very important here. 

• With increasing T, polarization tends to continuously penetrate into the superfluid core of 
the trap, until at, or even below Tc, the polarization is uniformly distributed throughout the 
cloud. 

• This superfluid Sarma or breached pair phase is limited in the amount of polarization it 
can accommodate, especially near unitarity. This applies to both traps and homogeneous 
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Figure 12. - Typical behavior of Tc as a function of l/kpa for spin polarizations p — 0.05, 0.1, 0.16, and 
0.3. The Tc curve splits into two disconnected curves around p — 0.14. This figure should be compared 
with left panel of Fig.|5|for the unpolarized case. From Ref. [75]. 



sy terns. 

Figure[21presents a plot of Tc for a homogeneous system as a function of 1 /kpa for various 
polarizations. This figure should be compared with the left panel of Fig.|5] With this comparison, 
one sees at once that there can be no superfluidity in the deep BCS regime, once the polarization 
is different from zero. And when superfluidity first appears on the BCS side of resonance it is 
associated with two T^s at given l/kpa. The larger the polarization p, the harder it is for a 
homogeneous system to support superfluidity, except in the EEC regime. 

To understand the meaning of these two Tc's we plot the superfluid density ns{T) as a func- 
tion of T in Fig.[Oj for several different values of l/kpa. If one focuses on the unitary case, 
for definiteness, one can see that Ug vanishes at two different temperatures. The lower Tc cor- 
responds to the onset of superfluidity. At temperatures below this, the (breached pair or Sarma) 
state is unstable. Similarly, at the upper Tc, superfluidity disappears in the usual way; it is de- 
stroyed by thermal excitations. For positive values of l/kpa, the figure shows that the ns{T) 
curves stop abruptly. This is indicated by the dotted segments of the curves which represent 
thermodynamically unstable solutions. 

The fact that Tc may be double-valued could have been anticipated in a simpler set of cal- 
culations performed at the strict mean field level, and discussed elsewhere in the literature [74]. 
Here one solves Eqs.( l2.Hl and ( I2.12l i only, without imposing ( I2.20l i and M.21\ . In Fig.^^ we 
present a plot of T*^^ as a function of l/kpa for a range of p. This quantity can be viewed as the 
pairing onset temperature T*, when there is a stable superfluid phase. In the inset of Fig.[T4lwe 
plot A(T) at different l/kpa forp = 0.3. Forp < 0.9 and sufficiently low T/^^, we find that 
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there are two T*^^ values. This structure implies that A is nonmonotonic [74] in T, as indicated 
by the bottom curve in the inset of Fig. O The two zeroes of A represent the two values of 
T^^^ . In contrast to the more conventional behavior (shown in the top curve for stronger pairing 
interaction), A increases initially with T at low temperature when 1/kpa is sufficiently small. 
We summarize by noting that these results indicate that temperature enables pairing associ- 




\lkpa 

Figure 14. - Mean-field behavior of T^'^ as a function of 1/fcFffl for different spin polarizations p. Shown 
in the inset is the pairing gap A(r) at different l/kpa for p — 0.3 which can vanish at two distinct 
temperatures. Here Ef = ksTF = f^k%/2m is the noninteracting Fermi energy in the absence of 
polarization. The quantity T^'^ represents the pairing onset temperature T*. 
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Figure 15. - Spatial profiles of the Gap (A), order parameter (A^c), and the density of the up (ri|) and down 
{ni) spin components and their difference {5n) for a unitary gas in a trap at T /Tf ~ 0.1, 0.14, 0.26, and 
0.28, from left to right. Except at very low T, most of the polarization is carried by fermions co-existing 
with noncondensed pairs, i.e., in the pseudogap regime where Asc = 0, but A 7^ 0. Here Tc/Tf ~ 0.27. 



ated with the breached pair or Sarma state. This was also inferred in Refs. [49] and [74]. In 
general superfluids, one would argue that these two effects compete. We may view, then, this 
unusual polarized phase as an "intermediate temperature superfluid". 

We now turn to the behavior of these superfluids in a trap. The same calculations are applied 
to the Sarma or breached pair state using the LDA to incorporate trap effects. Figure^jshows the 
resulting behavior at unitarity for polarization of 15% and for various temperatures from below 
to just above Tc. The upper panels plot the order parameter Asc and the (total) gap parameter A. 
Superposed on these plots is the polarization 6n. The lower panels present the density profiles 
for each spin state. Several important features can be gleaned from the upper panels. At the 
lowest temperatures the bulk of the polarization is in a region where Asc = 0, but A 7^ 0; 
thus polarization is excluded at low T from the superfluid core. Moreover, it can be seen that 
an excitation gap A is present throughout most of the cloud. Whenever A 7^ one can infer 
that both spin states are present. For non-zero A, the particle profiles are necessarily different 
from those of a non-interacting gas. The bulk of the polarization appears in the "normal, mixed 
region", and within this portion of the trap there are strong interactions between the two spin 
states. Only at the very edge of the trap is there an exclusively majority component (and here 
A = 0). This region is occupied by a non-interacting Fermi gas and can thereby be used to set 
the temperature scale for the trapped cloud. 

We end this Section with some comparisons with experiment. We address qualitative effects 
in Fig.^]and semi-quantitative effects in Fig.^] The former is for the unitary case and the latter 
is in the BEC regime. The left panels in Fig.ll6lshow data from the Rice group Ref. [77]. The 
upper figure plots the density profiles for each spin state and the lower panel, their difference. 
The unpolarized core is evident, as is the sharp edge beyond which polarization abruptly appears. 
This behavior has been interpreted [21,77] as suggesting phase separation. 

Because the aspect ratios for the trap are not maintained in the profiles, it has been argued 
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Figure 16. - Comparison of tiieoretically computed density distributions of the up and down spin compo- 
nents (upper row) and their difference (lower row) at unitarity with the Rice group data. The central issue 
here is not the quantitative comparison, but rather that (as seen in the theoretical curves) polarization is 
confined to outside the superfluid core. In the Rice data, has a narrower distribution due to non-LDA 
effects. Data from Ref. [77]. 



[21,77] that the LDA scheme may not be appropriate for addressing experiments on these highly 
anisotropic traps. In the experimental data, non-LDA effects lead to a narrower distribution for 
n| along the axial direction while it is broadened in the radial direction. Despite this caveat we 
plot our counterpart theoretical results at T = 0.06Tp (right panels) for qualitative comparison. 
This plot is designed primarily to introduce theoretical observations (which can be superposed, in 
effect, on the experimental plot) concerning where the superfluid phase resides in the trap. This 
knowledge cannot be directly gleaned from these particular experiments. These low T theoretical 
results are rather striking, for they make it clear that when the order parameter is present the 
polarization is excluded. Rather, the polarization appears in the Fermi gas region outside the 
condensate core. At higher T and lowerp (e.g, T — O.lTp, p = 0.15), the polarization is carried 
largely within the "pseudogap" regime, where there are strong pairing correlations, but no long 
range order. 

We turn finally to a comparison of data from the IVIIT group [22], and in the BEC regime. The 
experimental data (for polarizations near 0.6) are plotted on the two right panels which show the 
particle density profiles for each spin state and their difference 6n. The upper panels correspond 
to temperatures which are believed to be in the normal state, whereas, the lower panels are for 
the superfluid phase. The lightly dotted curves in the experiment represent extrapolations of 
the Thomas-Fermi fits to the curves at large radii. In the superfluid phase, the data show that the 
polarization at the core center is considerably smaller than it is for the normal state. One apparent 
difference between above and below Tc experiments, is that at low T (i) the minority component 
has contracted into the center of the trap, (ii) Another signature is a kink in the majority profile. 
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Figure 17. - Comparison of theoretically computed density distributions of the up and down spin compo- 
nents and their difference Sn at dilferent temperatures in the BEC regime with MIT data from Ref. [22], 
for p — 0.6. Upper/lower panels are above/below Tc. The arrow in the bottom theoretical curve for Sn 
indicates the condensate edge which is reflected in features in both and n^, much as in the data. 



Yet another signature is that (ill) there is a clear bi-modahty in the minority component. 

This behavior can be compared with the theoretical results plotted in the left two panels 
above and below Tc for a polarization p = 0.6. The kink in the majority component at low T 
can be clearly attributed to the edge of the condensate as denoted by the vertical arrow (where 
= 0). The bi-modality in the minority component is amplified in the lower inset to make it 
more evident. These theoretical plots, at a qualitative level, exhibit the three central features of 
the data noted above. 



7. - Conclusions 

Throughout this Review we have stressed that temperature is important in understanding and 
characterizing the ultracold Fermi gases. Experimentally, one is never strictly at T = 0; as a 
result there are thermal excitations of the gas which need to be characterized both theoretically 
and experimentally. At a formal level, we have discussed how temperature enters into the stan- 
dard BCS-BEC crossover theory by leading to a distinction between the superconducting order 
parameter A^c and the gap parameter A. This difference reflects the existence of pre-formed 
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pairs above Tc and noncondensed pairs below. 

At a more physical level we have shown how temperature changes the character of the gas 
from a superfluid at low T to an unusual (but strongly interacting) normal fluid at moderate T, 
and to ultimately a free Fermi gas at high T. This strongly interacting normal fluid is most inter- 
esting, for unlike the Landau Fermi liquid (or gas) which exists just above Tc in a strictly BCS 
superfluid, here there is a normal state excitation gap or pseudogap. By looking at three different 
experiments, we have provided evidence for this normal state excitation gap (in RF pairing gap 
spectroscopy [62] and in thermodynamics [8] and for the below-Tc counterpart (noncondensed 
pairs), via the shape of the particle density profiles [68, 69]. 

In this Review we have discussed how temperature can be measured through adiabatic sweep 
thermometry [3] and through Thomas-Fermi fits to density profiles [68, 8]. We have shown how 
temperature played an important role in establishing the first generation and earliest evidence 
for superfluidity - based on fast sweep experiments which yield [3] a phase diagram in the 
temperature vs 1/ /c^a plane. It also played an important role in the second generation indications 
for superfluidity associated with thermodynamic al measurements [8]. 

Finally, we have addressed a new and exciting class of experiments [20-22] involving spin 
population imbalance and shown that here too temperature plays a critical role. From an exper- 
imental point of view, the most exciting features of these experiments are that they show how 
the density profiles can be used to establish the transition temperature Tc- Moreover, through the 
wings of the profiles, they provide a theory-independent mechanism for thermometry. From our 
theoretical perspective, which focuses on temperature [75,56], what is also exciting is that (i) the 
stability and character of these polarized superfluids is very sensitive to temperature and (ii) even 
in the normal state one sees strong interactions between the two spin components, which we as- 
sociate with the finite T pseudogap effects we have been discussing in this Review. As stated in 
Ref. [22], "Akeady at high temperatures, above the phase transition, the larger [majority] cloud's 
profile is strongly deformed in the presence of the smaller [minority] cloud, a direct signature of 
interaction." 

Our discussion, thus far, has not addressed phase separation [42] or the exotic LOFF states 
[23, 25] associated with polarized superfluids. The latter, which correspond to a condensate 
with a finite momentum q, may well be ground states of these spin polarized gases for some 
range of scattering lengths a, since the polarized BCS like phase (with q = 0) is only stable at 
intermediate temperatures. These LOFF phases are discussed elsewhere in this volume. There 
are arguments in the literature that in some more general form they play an important role near 
unitarity [48, 40] in trapped gases. At a theoretical level, LOFF states may be included [55] 
following the formalism we have outlined in this Review. 

Underlying the interest in this general field of ultracold gases is the possibility that they may 
shed light on the high temperature cuprate superconductors either in the context [17, 19] of the 
BCS-BEC crossover scenario or in the context of optical lattices and Hubbard model simulations. 
For the former, one can return to the question of where they would fit on a phase diagram plot 
(such as that presented in Fig. |5] which addresses Tc vs l/kpo,)- When the calculations are 
properly redone for c?-wave pairing on nearly two dimensional lattices, the values and shape of 
the Tc curve are in quite good agreement with experiment [16-18]. Moreover, Tc vanishes well 
before the BEC limit is reached [26]; one can, then, deduce that in this scenario the cuprates are 
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close to the unitary regime, just as are the ultracold gases. 

Whether or not this crossover picture turns out to be correct for high Tc, another very im- 
portant rationale for its study is the possibility of generalizing what is undoubtedly the most 
successful theory in condensed matter physics: BCS theory. We now understand from the ultra- 
cold gases that the nature-made superconductors to which the original theory has been applied, 
are only a very special case of a much more general class of superfluids. And it is most fitting 
on this fiftieth anniversary of BCS that we pay homage to this most remarkable of theories by 
recognizing its even greater generality. 

* * * 
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